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Abstract 

This paper calculates the exact quadratic variation in space and quartic variation 
in time for the solutions to a one dimensional stochastic heat equation driven by a 
multiplicative space-time white noise. 
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1 Introduction 



We consider solutions (X t (x) : t > 0, x G R) to the one-dimensional stochastic 
heat equation 

dX = a AX dt + a{X) dr] (1) 

where a > 0, a : R — > R and rj(dx, dt) is a space-time white noise on 
R x [0,oo). Such equations arise in several settings, for example generalised 
Edwards- Wilkinson models for the roughening of surfaces [1], continuum lim- 
its of particle processes (see for example [2], [4], [8]), and continuous space 
parabolic Anderson models [3]. 
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Under reasonable conditions on a, there exist continuous solutions which are 
known to be Holder continuous of order | — e in space and of order | — e in 
time, for any e > 0. The aim of this note is to calculate the exact variations. 
We show that the paths x — > X t {x) have finite quadratic variation, for fixed 
t > and over an interval [Ai, A2], given by 

Jim £ (X t ( Xj ) - X t (x^)f = i- J2 a 2 (X t (x))dx, 

and the paths t — > X t {x) have quartic variation, over an interval [Ti, T 2 ], given 
by 

Here {xj} and {t,} are equally spaced partitions and the convergence is in 
probability. Precise statements are given later in Theorems 4 and 6. The results 
are local in character, being dependent only on the nature of the noise and 
the small time asymptotics of the Green's kernel, and we will remark later on 
similar results that hold for related stochastic PDEs. In the case of constant 
a the variations can be checked by an L 2 argument using the Gaussianity of 
the solution. For general a the proof involves showing that the increments are 
approximately conditionally Gaussian. 

The result is interesting when considered as a method for parameter estima- 
tion since the value of the drift parameter a can theoretically be found from 
an arbitrarily small segment of a single path. For stochastic ODEs the drift 
parameters are typically not reflected in the local structure of paths, since 
changes in value of the drift produce absolutely continuous changes in the law 
of the process on finite time intervals. Here however, exploiting say the time 
variation, we may choose 

, 3(T 2 -T 2 ) TU°\X tj {x)) 

Oi — - 

™ £J y(x tj {x)-X t M)Y' 

as an estimator of a, and an analogous estimator is possible that exploits the 
spatial variation. Consideration of the estimates involved in our results sug- 
gests, at least when a is smooth, bounded and bounded away from zero, that 
i?[|Qi n — ot\ A 1] = 0(r?T 3 / 20 ). Details on this and other parameter estimation 
problems will be found in the forthcoming thesis [6]. 



2 The linear case 

In this section we suppose that a(x) is a constant. By a linear change of 
variable we may suppose that a(x) — a — 1. The solution of (1) is then 
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defined as 

XAx) — \ G t (x — u)X (u)du + / / G t - r (z — u) r](du, dr) (2) 
Jb, Jo Jr. 

where G t (x) = exp(— x 2 /At) is the heat kernel. If we suppose some growth 
rate on the initial condition X , then the first integral on the right hand side of 
(2) is a smooth function of x G R and t > and will not affect the variations 
we consider. Without loss of generality we therefore assume for the rest of this 
section that X (x) = 0. 

The covariance structure of the solution (2), with X = 0, can be calculated, 
using the isometry for space-time white noise integrals (see Walsh [9] chapter 
2) and the semigroup property of the heat kernel, as, for s < t, 

E [X t (x)X s (y)] = [ r= 1 exp (- J*~ y) * ) dr. (3) 



2.1 Spatial variation 



Throughout this subsection we fix t > and A\ < A2. We consider the 
quadratic variation of x — > X t (x) over the interval [A±, A2]. From the covari- 
ance structure (3) we may deduce, by integration by parts, 

E[X t (y)X t (x)} 

FT ( (y-x) 2 \ (y-x) fy-x\ \y-x\ 



where erf denotes the error function defined by erf(x) = f£e * 5 ' dt. In 



particular, taking x = y, we have E[X 2 (x)] = \jt/2n. Define, for 5 > 0, the 
jointly Gaussian increments 

A(x, 5) =X t (x + S)- X t (x), A(y, 5) = X t (y + 5)- X t {y). 

From the covariance structure (4) we have that 

* /f exp (Jiy- e rf (^y-{ + Oify (5, 

Here we have expanded in small 5 and used the expansion erf (5) = -^5+0 (5 s ). 
Furthermore, supposing that x + 5 < y, we have 
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E[A(x,5)A(y,5)} 

= E [X t (x + 5)X t (y + 5)-X t (x + 5)X t (y) - X t (x)X t (y + 5)+ X t {x)X t {y)\ 




2 exp 



(y - xf 



— exp 



St 
St 



exp 



+ 2^erf 



(y-x + 5) 2 
St 
y-x 



2V2t 



exp 



(y-x + S) erf / y - x + 5 \ _ (y - x - 5) [y - x- 6" 



(y ~ x? 



4 V 2V2t 

Using the Taylor expansions 



St 



exp 



St 



2V2t 



At 



erf ( V^±l ) = erf ( V^£\ ± _^ exp 



(y - ^ 

8t 



+ 0(5 2 ). 



2V2t J \2V2lj V2nt 

we conclude that 

E[A(x,6)A(y,6)}=0(5 2 ). 

Here we have expanded in small 5 and used t,t~ x ,x, y = 0(1) where appro- 
priate. The expression 0(5 2 ) denotes a quantity bounded by C5 2 , for small 
5, with the same constant C for all x,y £ [A 1 ,A 2 ] satisfying x + 5 < y. We 
will need two higher Gaussian moments which can be read off from the lower 
moments (for example using E[X 2 Y 2 ] = E[X 2 ] E[Y 2 ] + 2(E[XY]) 2 ). When 
x + S < y we have 



E 



A\x,5)}=0(5 2 ), 
5 2 



E 

A 2 (x,5)A 2 (y,S) 



+ 0(5 6 ). 



(6) 
(7) 



Proposition 1 Define, for j = 0,1, ... ,n, a space grid by Xj — A ± + j5, 

where 5 = \{A 2 — Ai). Then the following limit holds in mean square: 

E (Uxj) - X t {x^)) 2 = \{A 2 - A,). 



PROOF. The desired L 2 convergence can be rewritten as 



hm y y E 



j=l k=l 



'a 2 (x 3 ,5)- 5 -) U 2 (x k ,5)- S - 



= 0. 
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Consider first the diagonal terms when j = k. These give a contribution 
bounded by 



A 4 (x v 5)+5A 2 (x j ,5) + 



S 2 



< nO(5 2 ) ^ as n — > oo 



using the estimates from (5,6) and the fact that 5 = 0(l/n). The off-diagonal 
terms give 



2EE £ 

j=i k=j+i 



A 2 ( Xj ,S)A 2 (x k ,S) - 5 -(A 2 (x v 5) + A 2 (x k ,5)) + j 



= 2E E 0(5 3 ). 
i=i fc=j+i 



The cancellation in the last equality uses the expansions in (5,7). The final 
sum converges to zero as n — > oo completing the proof. □ 



2.2 Temporal variation 



Throughout this subsection we fix x G R and times < 7\ < T 2 . We con- 
sider the quartic variation of t — > X t (x) over the interval [Ti,T 2 ]. From the 
covariance structure (3) we deduce, for s <t, 

E [X t (x)X s (x)} = (VtT~s - Vt^) . 

Suppose that s,t 6 pi,T 2 ] and < 5 < t — s. Define the jointly Gaussian 
increments 

A(s, 5) = X s+S (x) - X s (x), A(t, 5) = X t+S (x) - X t (x). 
Using the covariance structure above we find 



E 



A 2 (s,5) 



1 / 



2^ V 



( ^2(s + 5) + \/Zs - 2V2s + 5 + 2^6) = —j= + 0(5 2 ). 



Furthermore 

E[A(s,5)A(t,5)] = 



2^ 



2y/t + s + 5 - Vt + s + 25 - Vt + s 



+ 2y/t - s - y/t-s + 5 - y/t-s-5 

= 0(5 2 ) + 0((t- s y 3/2 5 2 ). 

Here we have expanded as a series in S and (t — s) _1 using that s, t, s+t = 0(1) 
where appropriate. We deduce from this covariance structure the following 
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higher asymptotics, which hold uniformly over s,t G [Ti,T 2 ] satisfying 5 < 
t — s, 



E 
E 



35 



A 4 (s,5) =- + 0(5 5 / 2 ), 



71 



E 



A s (s,5)}=0(5 2 ), 
9^ 

7T 2 



+ 0(5 7/2 ) + 0((t- s)- 3 5 5 ). 



(8) 
(9) 
(10) 



We have used the expansion of Gaussian higher moments in terms of the 
covariance function, for example E[X A Y A ] = 9(E[X 2 ]E[Y 2 ]) 2 + 2A(E[XY]) A + 
72i?[X 2 ]i?[y 2 ](i?[Xy]) 2 , keeping only the leading terms of the asymptotics in 
5 and (t — s) _1 that we will need. 

Proposition 2 Define, for j — 0, 1, . . . , n, a time grid by tj = T\ +j5, where 
5 = ^(T 2 — Ti). Then the following limit holds in mean square: 



n 4 3 

E - X^xj) = -(T 2 - Ti). 



PROOF. The desired L 2 convergence can be rewritten as 



hm y y E 



i=l fc=l 



0. 



Consider first the diagonal terms when j = k. These give a contribution 
bounded by 



HE 



A%,5) + 6 ^A%,5) + 9 4 

71 7T Z 



< nO(5 2 ) ^ asn^oo 



using the estimates from (8,9) and the fact that 5 = 0(l/n). The off-diagonal 
terms give 



n n 



2E E e 

j=l k=j+l 



A%, 5)A A (t k , 5) - -(A 4 (t,, 5) + A A (t k , 5)) + % 

71 7T Z 



2E E {(k-j)- s O(5 2 ) + 
j=l k=j+l 



The cancellation of the moments in the last equality uses the expansions from 
(8,10). The final sum converges to zero as n — > oo completing the proof. □ 
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3 The non-linear case 



We consider the equation (1) with a globally Lipschitz and satisfying a linear 
growth condition. By linear time scaling we may assume again that a — 1. We 
suppose that X is in the following space C tem of functions with slower than 
exponential growth 

C tem = {/ : R - R : ||/|| A := sup \f(x)\e~ Xx < oo, for all A > 0}. 

X 

Give Ctem the topology induced by the seminorms ||/||a for A > 0. The choice of 
C tem is not crucial but allows us to quote the following results. Let rj(dx, dt) be 
a space-time white noise on R x [0, oo), adapted to a filtered probability space 
(fi,^ 7 , {T)t,P)- In Shiga [7] Theorem 2.1 it is shown there exists a pathwise 
unique adapted, continuous C tem valued solution, which moreover satsifies the 
moment bounds, for any p, A, T > 0, 

E[\X t (x)\ p ] < C(p,\,T)e x W for all x e R and * e [0,T], (11) 

and the increment estimates 

E [\X t (x) - X s (y)\*>] < C(p, A, T)(\x- yf 2 + \t - sf/ 4 ) e^M+M) (12) 

for all x, y e R and s,t G [0,T]. In fact these, more or less standard, bounds 
are not explicitly stated in [7] but are implicit in the proofs found in the 
appendix. 

The solution can be written as X t (x) = J G t (x — u)X (u)du + X t (x) where 
X t (x)= f I G t - r (x -u)a(X r (u))r)(du,dr). (13) 

JO JR 

Since (t, x) — > / G t (x — u)X (u)du is smooth on t > 0, the quadratic and quar- 
tic variations of X t (x) will exactly coincide with those of X t (x). We therefore 
now restrict our consideration to the process X t (x). We will calculate the spa- 
tial and temporal variation in subsections 3.1 and 3.2 respectively. In 3.3 we 
describe some methods to transfer these results to related equations. 

3.1 Spatial variation 

We fix t > and consider the spatial increments defined by 

A(X,x,5)=Xt(x + 5)-X t (x). 

The key idea is that the main contribution to A(X,x,5) comes from the 
noise near time t. We will approximate the increment A(X,x,S) by the term 
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a(X t{s) (x))A(x, 5) where t(5) = t - 5 4 / 3 and 



A(x,5) = / / (G t - r (x + 8 — u) — G t - r (x — u)) r](du, dr) 

Jt(S) J 

ft(S) f 

+ J J {G t - r { x + 8 — u ) — G t - r { x — u)) r](du,dr) 

for fj(dx, dt) an independent space-time white noise. The integral against 
fj(dx,dt) is included so that A(x, 5) has the same Gaussian distribution as 
A(x,S) from the linear theory. Moreover A(x, S) is independent of u{X s : s < 
t(5)}. The value of t(5) is chosen to optimise the estimate in the following 
lemma, which shows the above approximation is valid. 

Lemma 3 For any A,T > there exists C(A,T) < oo so that 



E 



A(X,x,5)-a(X m (x))A(x,5) 



< C(A,T) <5 4/3 e AN 



for allx G R ; < 5 4 ^ < t A 1 andt<T. 



PROOF. We split the difference into three parts: 

A(X,x,S)-a(X m (x))A(x,5) 
t 

= J J (G t - r (x + 5 - u) - G t - r (x - u)) (a(X r (u)) - a(X t(5 )(x))^ r](du,dr) 



t(5) 

t(8) 

+ J J (G t ^ r (x + 5 — u) — G t - r (x — u)) a(X r (u)) r](du, dr) 



t(S) 

- J J (G t - r (x + 5 - u) - G t - r (x - u))a(X t( s)(x))fj(du,dr) 



= 1 + 11- III. 

We estimate each of these terms separately. We need the deterministic esti- 
mates, for A > 0, 



J* j G 2 t „ r {x-u)e XM dudr < C(A, T)(t - s^V 1 * 1 



(14) 



and 



J S J (G t _ r (x + S-u)- G t _ r {x - u)f e m du dr 

<C(\,T)5 2 (t-s)- 1/2 e x W, (15) 
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valid whenever 5 G [0, 1], < s < t < T and x e R. These (and the similar 
estimate (20) below) follow from straightforward calculations and we omit the 
proofs. Using the moments (11) and estimate (15) we have 









E 


\Ilf 


-ITS 



{G t _ r {x + 5-u)- G t _ r {x - u)f E a 2 (X r (u)) 



du dr 



< C(X, T)e AN <5 2 (t - t{5)y 1/2 = C(X, T)e AN 5 4/3 . 

The estimate on E^iTi"! 2 ] is entirely similar. Using the increment bounds (12) 
we bound -E^/l 2 ] by 

t 

C(A, T)e AN / / {G t - r {x + 6-u)- G t - r (x - u)f (\x -u\ + 5 2/3 ) e XH du dr 



t(S) 



< t _ r {x — u) [\x — u\ + 

rt 
lt(5) 



5 + 5 2/3 ) e XM dudr 



< C(A,T)e 2A|x| <5 4/3 + C(A,T)e A|:c| f f G 2 t _ r {x - u)\x - u\e XM dudr 

< C(X, T)e 2AN 5 4/3 + C(X, T)e AN /' / G t - r (x - u)e XH du dr 
<C(X,T)e 2X ^ 3 . 

The second inequality here uses (14) and in the third we used the bound 
sup tz \zG t (z)\ < oo. Noting A > is arbitrary, we see that combining the 
bounds on the terms J, 77 and III completes the proof. □ 



Using this approximation we can establish the spatial variation in the non- 
linear case along the same lines as for the linear case. 

Theorem 4 Fix A\ < A 2 and t > 0. Define, for j — 0, 1, . . . , n, a space grid 
by Xj — Ai + j5, where 5 = \{A 2 — Ai). Then the following limit holds in 
probability: 

n ~ 2 1 
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PROOF. We break the required convergence into three parts as follows. 

1 f A 2 



£A 2 pr,^,5) / a 2 (X t (u))du 

j=l * Ja i 

J2 (A 2 (X, Xj ,5) - a 2 (X m ( Xj ))A 2 (x j: S 



+ ^AX m (x J ))[A 2 (x J ,5) 



2 / 

1 rM 



(16) 
(17) 



+ 2 EX P^)) - 2 ^W"))**- (18) 

The third term (18) converges almost surely to zero as n — » oo using the 
uniform continuity of (s,x) — > X s (x) and a Riemann sum approximation to 
the integral. The first term (16) converges to zero in L*(P) since 

± E [| A 2 (X, Xj , 5) - a 2 (X m (x 3 ))A 2 (x 3 , 5)\ 
j'=i 

< £ [|A(X, *,-,*) - a(X t(5) ( a ;,))A(^,5)| 2 ]) 1/2 
j'=i 

• (E [|A(X,a;,,5) + ( T(X^ ) (^))A(a; J , ( 5)| 2 ]) 1/2 

< C(t,A 1 ,A 2 )5 2 / 3 J2{ E [A 2 (X,x v 6) + a 2 (X m (x J ))A 2 (x J ,5)]) 1/2 

< C(t,Ai,A 2 )(J 2 / 3 nO(<J 1 / 2 ) ->0. 

The first inequality is Cauchy-Schwartz, the second inequality uses the bound 
from Lemma 3 and in the final inequality we have applied the moments from 
(11) and (12) and the independence between A(xj,S) and a{X t {S){xj)). 

We will now show the second term (17) converges to zero in L 2 (P) by mim- 
icking the linear case. Indeed, using the independence between A(xj,5) and 
cr(X t (8)(xj)) and the moments (5,6), we have 



E 



±a\X m {x 3 ))[A 2 {x^)- 5 - 



< E 



^\X m {x 3 ))0{5 2 ) 

3=1 



which converges to zero using the moments bounds (11). Similarly, using (7), 



E 



± a 2 {X t{s) { Xj ))a 2 {X m {x k )) (a 2 (x„6) - 5 -) (A 2 (x k ,5) - f 



= E 



E ^{X^x^X^x^W) 
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which also converges to zero. This shows the second term (17) converges to 
zero in L 2 (P) and completes the proof. □ 



3.2 Temporal variation 



Fix x and define temporal increments via 

A(X,t,5)=X t+s (x)-X t (x). 

The key idea is that the main contribution to A(X, t, 5) comes from the noise 
near time t. We will approximate the increment by cr(X t (s)(x))A(t,5) where 
t(5) = t - <5 4 / 5 and 



A(t,5)= [ [ (G t+ s- r (x - u) - G t - r (x - u)) r)(du, dr) 

Jt(S) J 

rt(S) r 

+ J j (G t+ s-r(x -u) -G t - r (x -u))fj(du,dr). (19) 



where fj(dx, dt) is an independent space-time white noise. We use here, and 
below, the convention that G t (x) = for t < 0. Note again that A(t,S) 
has the same Gaussian distribution as A(t,8) from the linear theory and is 
independent of a{X s : s < t(5)}. The value of t(5) is chosen to optimise 
the estimate in the following lemma, which shows the above approximation is 
valid. 



Lemma 5 For any A, T > there exists C(X,T) < oo so that 



E 



A(X,t,5)-a(X t{s) (x))A(t,5) 



< C(\,T)5 4/5 e 



an 



for all x e R, < 5^ 5 < t A 1 andt<T. 
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PROOF. We split the difference into three parts: 

A(X,t,5)-a(X t{s) (x))A(t,5) 
t+s 

= J J (G t+ s-r(x - u) - G t -r(x - u)) (a(X r (u)) - a(X m (x))^j r}(du, dr) 



t(S) 



+ J J (G t+ 5- r (x - u) - G t - r (x - u)) a(X r (u)) r)(du, dr) 



t{5) 

- J J (G t +s-r(x - u) - G t _ r (x - u)) a(X t{5) (x)) fj(du, dr) 



= 1 + 11 - III. 

We estimate each of these terms separately. We need the deterministic esti- 
mate, for A > 0, 

f J (G t+ s-r(x - u) - G t ^ r (x - u)f e x ^du dr 

<C(X,T)5 2 (t-s)- 3 / 2 e x ^, (20) 

valid whenever S e [0, 1], < s < t < T and x G R. Using the moments (11) 
and estimate (20) we have 



E 



\II\ 



rt(S) r „ r 

J J (G t+S - r (x u) G t _ r (x - u)f E [a 2 (X r (u)) 



du dr 



< C(X, T)e x ^5 2 (t - t{5)y 3/2 = C(A, T)e x ^5 4/5 . 



The estimate on i?[|TZT| 2 ] is entirely similar. Using the increment bounds (12) 
we bound E [\I\ 2 ] by 



C(\,T)e x ^ J J (G t+ s- r (x-u)-G t - r (x-u)) 2 (\ x — u\ + du dr 



which, as arguing as in Lemma 3, is bounded by C(A, T)e 2A ' x '<5 4//5 . Combining 
the bounds on the terms J, II and III completes the proof. □ 



Using this approximation we can establish the spatial variation in the non- 
linear case along the same lines as for the linear case. 

Theorem 6 Fix < Ti < T 2 and x G R. Define, for j — 0, 1, . . . , n, a time 
grid by tj = Ti + j5, where 5 = ^(T 2 — Ti). Then the following limit holds in 
probability: 

J™ t - ^-MT = I C Ax t (x))dt. 

j=l 1 



12 



PROOF. We use one more trick in the approximation for A(X,tj,5). We 
approximate this increment by a(X t .^(x))Aj(tj, x) where tj(5) = tj — 5 A ^ 
and where Aj(tj, x) are defined as in (19) using an I.I.D. sequence (fjj(dx, dt) : 
j — 0,1,. . .) of independent space-time white noises. 

We break the required convergence into three parts as follows. 

3 r T 2 



" 1 fJ-2 

J2A 4 (X, tj ,5)-- <j\X t {x))dt 

= ± (A\X,t„5) - a^X^x^A^S 

3=1 

n ( 35\ 

3 r T 2 



Q X n r l 2 

71 p{ 7T 7Ti 



(21) 
(22) 
(23) 



The third term (23) converges almost surely to zero as n — > oo using the 
uniform continuity of (s,x) — > X s (x) and a Riemann sum approximation to 
the integral. The first term (21) converges to zero in L(P) 1 since 



E E [| A4 (^ $ - ^(Xw) {x))K%t 3 , 5) 
j'=i 

< CE (E [|A(X,t,,5) - ( x(X, 3( , ) ( a ;))A,(t,,5)| 2 ]) 1/2 

• (E [A 6 pr^,5) + ( T 6 (X t(5) (t,))A5(t J ,5)]) 1/2 

< C(|a;|,T 2 )(5 2/5 nO(5 3/4 ) ^0. 

The final inequality uses the bound from Lemma 3, the increment moments 
from (11) and (12) and the independence between A(tj,5) and a(X t .( S )(x)). 

We will now show the second term (22) converges to zero in L 2 (P) by mim- 
icking the linear case. We consider the double sum 



E^ 

j,k=i 



35' 



35 



a\X t]{5) {x))a\X tk{5) {x)) A%,5) - - Ai(t k ,5) - - 



71 



71 



. (24) 



For a term corresponding to j, k satisfying \j ; — k\ < n 1 / 2 we apply Holder and 
use the independence between A(tj,5) and a(X t .^)( x )) an d the asymptotics 
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(8,9) for moments of A(t, 5) from the linear case to find 
E 



a\X t](s) {x))a\X tk(s) {x)) Ufa, 8) - ^ j Ufa,*) - ^ 



< I E 

E 

0(5") 



o*(X th{S) (x)) Ui(t k ,S)-™\ 



1/2 



1/2 



This implies that the sum over j, k satisfying \j — k\ < n 1 ^ 2 in (24) converges 
to zero as n — > oo. Let G(k,j) = a{X s (x) : s < t k (5),x G R} V cr{fjj(dx, dt) : 
t > 0,x E R}. For k > j + n 1 / 2 , the variable A k (t k ,5) is independent of Q{k,j) 
whereas Aj(tj : 5) is G(k,j) measurable. For such j, k we therefore have 



E 

< E 



a\X tj{S) (x))a\X tk{S) (x)) Ufa, 8) - ^) Ui(t k ,8) 36 



TT 



35' 



a\X m (x))a\X thW (x)) A%,5) - - 0(5 5 / 2 ) 



TT 



= 0{5 9 ' 2 ). 



1/2 



For the first inequality here we conditioned on G(k,j) and used the asymp- 
totics from the linear case to calculate E[Aj,(t k , x) — (35 /it)}. This implies that 
the sum over j, k satisfying k > j + n 1 ^ 2 in (24) converges to zero as n — > oo. 
A similar estimate holds for the terms j > k + n 1 / 2 and this completes the 
proof. □ 



3.3 Related results 



The results above rely only on the correlation of the space-time white noise 
and the small t asymptotics of the heat kernel G t (x), and so hold for many 
variations of the stochastic PDE considered above. We comment here on some 
possible changes. 

The addition of a drift term such as 

dX = a AX dt + b(X) dt + a(X) d7] (25) 
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will not affect the results. Under Lipschitz and growth conditions on b the 
moments and increment bounds (11, 12) still hold true. In particular these 
imply that the solution (t,x) — > X t (x) is Holder continuous. The extra term 
in the solution J * / G t - r { x — u)b(X r {u))dudr is then continuously differentiable 
when t > 0, x G R, and therefore will not affect the variations studied here. 
This is also clear when the addition of the drift induces an absolutely contin- 
uous change in the law of the solutions (for example when b is bounded and 
cr(x) > <j > 0). Furthermore the results should go over for the equation (25) 
where b(X t (x)) and a(X t (x)) are replaced by adapted fields b(t,x) and a(t,x) 
provided that they are known to be Holder continuous in (t, x) and satisfy 
reasonable moment conditions. 

An example that does not fit into our Lipschitz assumptions above is the 
equation (1) with a(x) = x 13 , for (5 G (0,1), and where one considers non- 
negative solutions. In particular the case (3=1/2 occurs as the density of 
Dawson- Watanabe processes (see Konno and Shiga [4]). For these equations 
the moment and increment bounds (11,12) still hold, since the proofs rely only 
on the linear growth bounds on a. The other place the Lipschitz assumption is 
needed is in Lemmas 3 and 5 where it is used to get estimates on E[\a(X t (x)) — 
a(X s (y))\]. For this we may use the inequality 



\a 



f> - & /f < \ a _ 6 |2/9 < (1 _ 0)^/2 + p 5 -i(i-f))W\ a _ 6 f 



valid for a, b > and 7, 5 > (and derived from Holder's inequality). Together 
with the increment bounds (12) one reaches, for \t — s\ < <5 7 and A > 0, 



E 



\(X t (x) f - (X s (y)f\ 2 ] < C(/3, A, T) + \x - y\5'^^ 2 ) e x ^ + ^ 



for all x, y and s, t G [0, T\. Using this in the proofs of Lemmas 3 and 5, letting 
t — t(5) = 5 1 with 7 = 4/(2 + P) and 7 = 4/(4 + (5) respectively, the same 
variation results can be checked to hold. 

For the equation (1) over a finite interval x G [a, b] with, say, periodic or 
Dirichlet or Neumann boundary conditions, the same variation results should 
hold on (a, b) x (0, 00). This is intuitively clear since the small time behaviour 
of the Green's function on the bounded interval remains the same as that of 
the heat kernel on the whole of R. In the linear case, that is with constant 
coefficients, a simple way to deduce the results for bounded intervals from 
the whole space results is via a change of measure argument. In Mueller and 
Tribe [5] Corollary 4 there is an expression for the Radon Nikodym derivative 
of the law of X t {x) for a solution on (a, b) with periodic boundary conditions 
with respect to a solution on the whole space (with a periodically extended 
initial condition). Dropping to subsequence, if necessary, to obtain almost sure 
limits, the variation limits therefore still hold. (The proof that the expression 
is truly a change of measure is incorrect in [5], but a simple corrected proof is 
given, as an errata, in the same journal). 
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